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Cumulative observational evidence including recent findings confirm that the surface density of 
dark matter (DM) halos por> = To po where ro and po are the halo core radius and central 
density, respectively, is nearly constant and independent of galaxy luminosity for a high number 
of galactic systems (spirals, dwarf irregular and spheroidals, elliptics) spanning over 14 magnitudes 
in luminosity and of different Hubble types. We point out that this universal relation is of the 
same kind as the Larson relation for the average surface or column density pox> of interstellar 
molecular clouds of size ro (valid irrespective of the cloud type and composition over more than 
four orders of magnitude of scales 0.001 pc < ro < 100 pc). Remarkably, its numerical value 
fj,0D — 140 Mq/pc 2 = (18.6 Mev) 3 is approximately the same (up to a factor of two) in all these 
systems. A simple physical explanation of the Constance of pod based on the virial theorem and the 
extensivity property of the total energy is given. The action of self-gravity (both of baryonic and dark 
matter) must be the origin of the universal character of pod for all these cored bounded virialized 
systems, pod being a dynamical scale determined by gravitational clustering. We theoretically 
compute pod and the density profile for small scales from the linearized Boltzmann-Vlasov equation 
for self-gravitating DM. We solve this equation during the matter dominated era in the Gilbert form. 
We use the appropriate initial conditions including the primordial inflationary power spectrum and 
the evolution through the radiation dominated era. The theoretically derived density profiles pu n (r) 
turn to be cored and the resulting value for poo is remarkably close to the observations. Our 
density profile pu n {r) decreases as r _1_ns//2 for small scales where n s ~ 0.964 is the primordial 
spectral index. This scaling is in remarkably agreement with the empirical behaviour recently found 
observationally: r _1 ' 6±0 - 4 . The value of pod implies that the DM particle mass is in the keV scale. 
We consider both fermions and bosons as DM particles decoupling either ultrarelativistically or 
non-relativistically. Our results do not use any particle physics model and vary slightly with the 
statistics of the DM particle. This analysis provides further evidence for the mass scale of the (cold) 
dark matter particle in the keV range. 



I. OBSERVATIONAL EVIDENCES, THE VIRIAL THEOREM AND THE EXTENSIVITY OF ENERGY 

Among the scaling laws put forward by Larson [l| for interstellar molecular clouds in our Galaxy there is the relation 
between the cloud size ro and the average density po- 

ro Po = Pod 



ely a 

0.001 pc < ro < 100 pc. The values given in ref.pj are: 



where the column density pqd is approximately a constant over more than four orders of magnitude of scales 

in ref.pj are: 

Pod = 10.5 10 21 = 162 



pc 2 
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where m# 2 stands for the mass of the Hydrogen molecule, main constituent of the interstellar clouds. Recent data 
averaged over high density regions of Taurus give [2j 

Hod = 5.14 10 21 ^ = 80 

cm z pc z 

The mean density of structures in the ISM vary between 10 and 10 5 atoms/cm 3 , significantly above the mean ISM 
density which is about 0.1 atoms/cm 3 or 1.6 10 -25 g/cm 3 . 

Growing recent findings point towards the same relation for the dark matter (DM) surface density hod in galaxies 
[H EL 0- Namely, that the product pojj = tq po where r$ and po are the halo core radius and central density, 
respectively, is nearly constant, over a large number of galaxies of different kinds, tq varies by two orders of 
magnitude [J 0,0: 

0.3kpc < r < 30 kpc and 10~ 25 g/cm 3 < p < 6 x 10~ 23 g/cm 3 . (1.1) 

This finding relates to data sets (high quality rotation curves, kinematics, galaxy-galaxy weak lensing signals) for many 
galactic systems spanning over 14 magnitudes in luminosity and of different Hubble type, dwarf disk and sphcroidals, 
spirals, ellipticals. In spite of their different properties, hod in galaxies is essentially independent of their luminosity 
and mass. 

Remarkably enough, /i d takes in all these systems (molecular clouds and galaxy systems) approximatively the 
same value p], 0, 0, 0] 

Hod ~ 140 — | = 6400 MeV 3 = (18.6 Mev) 3 . (1.2) 

We give in this section a simple physical explanation of the Constance of hod and its value based on the virial theorem 
and the extensivity property of the total energy. We derive in the next section the value of hod from the resolution 
of the linearized Boltzmann-Vlasov equation for self-gravitating DM. 

Assuming that the matter distribution in the galaxy is characterized by a scale ro , the matter density can be written 

as 

p(r) = poF(j^j , F(0) = 1. (1.3) 

Algebraic fits to the cored density profile 00] and thermal profiles [l(| are particular examples of ea. (|1.3|) . We have 
for the Burkert [§], Spano and Sersic profiles [§] (denoted Fb, Fs p and Fs, respectively): 

1 1 _ — v 

Fb{x) = -— — r-— — 2\ , F Sp (x) = r , F s (x) =e xn , x = — , 

(l + .x) (l + x 2 ) p (l + x 2 )* r o 

where | < n < 10 fits most of the galaxies in the Sersic profile. The de Vaucouleurs profile is a particular case of the 
Sersic profile for n = 4. Notice that both the Burkert and the Spano profiles decay for large distances with the same 
1 /r 3 tail as the cuspy NFW profile [HJ 

A 

Fnfw(x) 



x (1 + x) 2 

The virial theorem for self-gravitating systems states that the total energy E is related to the average potential energy 
(U) and the average kinetic energy (K) by [HI 

E=\{U)=-{K) . (1.4) 
We can therefore express the total energy E in terms of the average gravitational potential energy as 

E = -\ G J |^ <p(r) p{r')) =-\G^rlJ |^ (F(x) F{ X ')) . 
Hence, the energy divided by the characteristic volume r-g goes as 

^ ~ G pi rg = G pl D . (1.5) 
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Extensivity of the energy requires E/r^ to stay fixed for large values of tq [13, HH, EH ■ Therefore, /iod must take the 
same constant value irrespective of the value of r . Eq. Q1.5P thus provides a simple explanation of the astronomical 
observations [3, 0, 0, 0] ■ 

We can estimate the average kinetic energy (K) as 

(K) = \ J d 3 r (p(r)) {v 2 ) = \ PQ rg (« 2 > | d 3 x ~ Po r 3 (« 2 > . (1.6) 

We thus find in general the average squared velocity (v 2 ) for any profile ea. (|1.3p from eqs. dl.4p . <| 1 . 5|) and (jl.6p . 

(« 2 ) ~ G mod r , (1.7) 

as in the particular case of the Burkert profile Q. Eq. (|1.7p is verified both for molecular clouds and galaxies. 

For the self-gravitating gas in thermal equilibrium the above argument eq. (|1.5p of the extensivity property of the 
energy leads to a density p(r) that scales as 1/r (lpj . 

The quantities tq and po depend on the particular galaxy (or molecular cloud) chosen and are therefore functions of 
the past history of the galaxy (or cloud). Instead, the product poo = ?*o Po given by eq. (|1.2[) is an universal number 
for all galactic systems and molecular clouds and hence independent of the previous history of the system, pqd can 
only depend on universal quantities. Since /ior> is the same (up to a factor two) for molecular clouds and galaxies, 
the action self-gravity (both of baryonic and dark matter) must be responsable of its value. As we show in the next 
section, /xod is a dynamical scale determined by gravitational clustering. We consider in what follows DM in galaxies. 



II. THE SURFACE DENSITY FROM THE LINEARIZED BOLTZMANN-VLASOV EQUATION 

The mass density P u n {r) can be evaluated theoretically solving the linearized Boltzmann-Vlasov equation for self- 
gravitating DM in the matter dominated (MD) era. It is convenient to recast such equation as a integral equation, 
namely the Gilbert equation which is a Volterra equation of second kind [l?], EH, Eil ■ To linear order in perturbations 
the distribution function of the decoupled particles can be written as 

f(x,p;t) = g /oO) + F 1 (x,p;t) 

where x,p are comoving coordinates, g is the number of internal degrees of freedom of the DM particle, typically 
1 < 9 < 4, fo(p) is the thermal equilibrium unperturbed distribution function at the decoupling temperature Td 
normalized by 



171 9 I 77^X3 hip) = Pdm = Pc , (2.1) 



d 3 p 

(27T) 5 

f2jvr = 0.233 is the DM fraction, p c is the critical density of the Universe 

p c = 3 M 2 Pl Hi = (2.518 meV) 4 , 1 meV = 1(T 3 eV (2.2) 

and Hq = 1.50 28 1 -42 GeV. Terms of order higher than one in iq are neglected in the Boltzmann-Vlasov equation. 
The physical initial conditions at t eq , the beginning of the MD era, are 

f(x,p;t eq ) =g /o(p)[l + S(x,t eg )] , Fi(x,p;t eq ) = g / Q (p) S(x,t eq ) , (2.3) 

where S(x, t eq ) are the density fluctuations by the beginning of the MD era. 
It is useful to integrate the fluctuations F±(p] x; t) over the momentum p, 

A(k,t)=m [ t^tt* F t (k,p-,t) where F l {k,p-t)=(d z xe- l3 ^F l {x,p;t). (2.4) 
J (27r) J J 

Its Fourier transform provides the matter density fluctuations in linear approximation today 

1 f 00 

Piin(r) = — — o — / kdk sin(fcr) A(fc,t t odav) , (2.5) 
Zir z r J 

where as customary we considered a spherical symmetric distribution. 
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We therefore have as initial conditions, 

A(Meg) = fi M Pc V 8{k,t eq ) , (2.6) 

where V ~ is the comoving horizon volume by equilibration. Namely, all fluctuations with k > k eq that were 

inside the horizon by equilibration are relevant here. More explicitly, k eq ~ 42.04 Hq = 9.88 Gpc -1 [25[ and 

1 f 

V ~ T3- ^ 773 where / ~ 1.346 10~ 5 . (2.7) 

The present linear treatment is valid for small scales k > k eq . Non-linear effects become important for large scales 
k < k eq and call for the use of the full (non-linear) Boltzmann-Vlasov equation or iV-body simulations. 

The density fluctuations 5(k,t) by the end of the RD era can be obtained analytically for subhorizon wavenumbers 
[25l l26l] . The initial conditions for S(k,t) in the RD era are the primordial inflationary fluctuations. With such initial 
conditions and solving the fluid equations for DM during the RD era, yields 25, 26] 



5(k,t eq ) = i A |</ia.| i 5 log 



Ay/2B k 
k eq (1 + V2) 2 



6 V2 - 15 \ = £ A \<t> k \ log ( 0.2637 B A ) , (2.8) 



keq 



where A ~ 9.6 and £? ~ 0.44 are constants that follow evolving the fluid equations 26], <fik are the primordial 
inflationary fluctuations of the newtonian potential [2(| l25j 

|0 fe |=V2^|A o | f-J , (2.9) 

|Ao| stands for the primordial power, n s is the spectral index, and ko the pivot wavenumber [20I [2l|. 

|A | ~ 4.94 10~ 5 , n s ~ 0.964 , fc = 2 Gpc" 1 . (2.10) 

It is convenient to define 



Then, the Gilbert equation takes the form [1 
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f u A(k u') 

A(k, u) — — / U[a (u — u')] — -7T du' — I [a u] where , 

a J [1-u'Y 

lf°° 1 f°° sin(yz) 

n[z] = — / dyyf (y) sin(yz) , I[z] = — / dy y f (y) and 

h Jo h Jo z 



p 2 k / 1 + z eq T d 

V — 77-T 1 z = au, a = — \ — — 

Td H V Sijwf m 



X 



1 



*2= / dyy z fo{y) , 1 + z e , = — ~ 3200 . (2.12) 



a 



<'ii 



u is a dimensionless time variable related to the scale factor by 

u = 1 - J — , a{u) = 0,69 , < u < utoday = 1 - v^eo — 0.982 , a(today) = 1 . 
V a (1 — u) z 

It follows from the resolution of the Gilbert equation eq. (|2.12[) that for late times the solution grows as [l9[ 

A(fc,i)^ day jjr(fc) (l + z eq ) (2.13) 

where T(fc) stands for the transfer function. A(fc, t) grows proportional to the scale factor in the linear approximation 
for all k < kf s - The free-streaming comoving wavenumber kf s increases with time as 1/\A + z. 
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Notice that T(0) = 1 and T(k — > oo) = 0. T(k) decreases with k according to the characteristic scale given by the 
free streaming wavenumber kf s where lf s — VE/kf s is the free streaming length '19]. The explicit expression of the 
comoving free streaming length is fl9j 

2^3 / 1 + Z eq 2 _ 1 2 



If. = °DM > °DM =3 < V '> . (2.14) 

c_dm stands for the primordial comoving squared velocity dispersion of the DM particles. That is, the velocity 
dispersion computed from the equilibrium thermal distribution fo(p) which can be expressed as 



(Tdm = \l TTT" — where h= I dy y 4 f (y) . (2.15) 



3 12 m 

It is convenient to introduce the variable 

7 = k run where r Hn = ^ = , (2.16) 

V3 fc /s 

and consider the transfer function T(k) as a function of 7. T(7) decreases by an amount of order one for 7 increasing 
by unit. Therefore, its Fourier transform pu n (r) eq. (|2.5[) . decreases with r having ru n as characteristic scale. 

T(k) is obtained by solving the Gilbert equation (|2.12ll [19] . We plot in fig. [T]T(fc) for Fermions (FD) and Bosons 
(BE) decoupling ultrarelativistically and for particles decoupling non-relativistically (Maxwell-Boltzmann statistics). 

The dark matter density eq. (|2.1|) can also be expressed as an integral over y as 

Pdm = ^gTlh. (2.17) 

The covariant decoupling temperature Td can be related to the effective number of UR degrees of freedom at decoupling 
gd and the photon temperature today T cm {, by using entropy conservation (28j : 

T d = 3 T cmb , where T cmb = 0.2348 meV . (2.18) 

We obtain the amplitude A(k,t) today by inserting eqs. (JUSJ, (pT7) . (JHHJ, ([2jl]) and (|2~T3]) into eq. ([2~TT|) for 

t — ^today 

A(fc, t today ) = ^ -^i n M f A (1 + z eq ) I A„| T(fc) ^-J log [c — J . (2.19) 

where c = 0.11604. Inserting eq. (|2.19p into cq. (|2.5[) yields the density profile today, 

Piin{r) = — ^ / |A | (fc r K „) -/ d7 7 s/ log sin 7 T( 7 ) , 

" ^ CT DM '"JO V ™eg r lin / \ ?"Zin / 

r K „ p K „(Q) = ^ ffi H ° / I A„| (fc r lin f- n ' /2 [°° *1 ^ lo S flT^M T ^ , (2-20) 

O 7T ® DM JO \ ^BQ '"Kw / 

It is very useful to relate the free streaming length to the phase-space density p/cr 3 [3, EH EH . p/c 3 is invariant under 
the cosmological expansion and decreases due to gravitational clustering. The phase-space density before structure 
formation (pdm / 1(7 dm) ano - today can be related as [lj| 

BL = \ E™ . p.^ 
CT s z "dm 



where the value today 



% ~ 5 x 10 3 = (o.is keV) 4 , (2.22) 

". (km/s) 



follows from dSphs observations [23|], and according to iV-body simulations 1 < Z < 1000 [24 
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Combining eqs. (|2.15p . (|2.17|) . (|2. 18[) and (|2.2ip we can express m and g& as 



^ Z ^n m = 0.2504 Z -Y % keV , (2.23) 



3 \/3 9 o-f A \ g 



L 2 1 2 



2 J 3 1 ^ 



S«* = -ZTT-^r !T7TfT2 -T (/2 ^ = 35 - 96 Z * 53 (/2 ^ ■ (2 - 24) 

For example, for fermions and bosons that decouple ultrarelativistically at thermal equilibrium eqs. (|2.23p and (|2.24|) 
yield [H 

/ Z\ 4 f 0.568 Fermions a 1 f 155 Fermions ,„ 

™ = W kcV \ 0.484 Bosons ' 5 ^ 54Z4 (l80 Bosons " ^ 

Notice that 1 < zi < 5.6 for 1< Z < 1000. 

We obtain the primordial DM dispersion velocity a dm from eqs. (|2.1[) . (|2.2p and (|2.2ip [l5l |. 

c dm = (3 ff 2 IW - ^) 3 (2.26) 

This expression is valid for any kind of DM particles. We find using eq. (|2.14h for the free streaming length, cq. (|2.26p 
for a dm j and e Q- (]2.22[) , 

It 269 5 /10\ 3 1 2 

r Hn = = -^r- kpc = 125.1 ( — ) kpc and — — = 2.34 10 13 Z? . (2.27) 

a DM 



V3 Z\ * \Z 



The velocity dispersion <tdm ~ 10 7 is very small since it does not take into account the self-gravity contrary to 
a s ~ 10~ 5 . (j dm is just the primordial velocity dispersion redshifted today. 
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The linearized Boltzmann-Vlasov equation provides a single galaxy configuration with characteristic size ru n ~ lf s . 
This is an universal length and it corresponds to the largest galaxy sizes eg . (| 1 . 1 [i . ru n > r . Non-linear effects would 
give a variety of galaxy configurations with smaller masses and sizes. 

Inserting ru n and 1 /o%m eq. I|2.27jl and n s eq. (|2.10p in eq. (|2.20|) yields for the density profile and the surface density: 

%n a /6 poo / \ 

p lin (r) = (5.826 Mev) 3 / log (eZ* 7 J sin I 7 — J T( 7 ) d 7 , (2.28) 

f°° / \ 

Mod = r K „ p im (0) = (5.826 Mev) 3 W 6 j 7 "=/ 2 log 7 J T( 7 ) d 7 . (2.29) 

where n s /2 - 1 = -0.518, n s /2 = 0.482, n s /6 = 0.160 and c = 43.6. 

We plot in fig. [2] pu n {r) / pu n (0) vs. X = r/ru n for Fermions (FD) and Bosons (BE) decoupling ultrarelativistically 
and for particles decoupling non-relativistically [Maxwcll-Boltzmann statistics (MB)]. We choose for the plots the 
value Z = 10 both suggested by TV-body simulations [2J] and by the observations for p, 0D [cq. (|1.2[) and Table I]. 
Anyway, the Z dependence of pu n (r)/ pn n (0) vs. x is very mild, through logZ, as shown by eqs. (|2.28[) - (|2.29[) . The 
displayed profiles are clearly cored, as expected, since T(k) decays for k > kf s ~ 1/nm- Moreover, the profile 
cq. (|2.28|) is flat at r = with a negative concavity at r = 0, namely p' Un (0) = and Pi' in {0) < 0. More explicity, 

n, (r) ^ t 2 Jn^T 2 ' 482 log(cZ3 7 ) T( 7 ) d 7 2 // /qn 

P' 1 "^^ 1 r< ^n™ i _ fL V I + o(x 4 ) = 1 + — '"^ ' + 0{x A ) 

PUn(0) 6 / °° 7 0.482 loghzi 7 ) T( 7 ) d 7 2 

We display in Table I the values of ru n pu n (0) and rf in PH„(0)/p;i,i(0) for the three particle statistics. We find that 
PUniS^BE > Piinifi)FD > Piin^MB- The more peaked density profile is the one for bosons and the more shallow is 
the non-relativistic one. The fermions profile being in-between the two other profiles. 




FIG. 2: The profiles pu„(r) / pn n (0) vs. x, where x = r/ru n for Fermions and Bosons decoupling ultrarelativistically and for 
particles decoupling non-relativistically (Maxwell-Boltzmann statistics). The bosons profile is the more peaked, the MB profile 
is the shallowest and the fermions profile is lying in-between. 



The value of the surface density r$ p(0) should be universal. This is why we identify ru n pu n (0) computed 
for a spherically symmetric solution of the linearized Boltzmann-Vlasov equation for self-gravitating DM with the 
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Particle Statistics 


Pod = rim pun(Q) 


rfin Plini ) 1 Plin(0) 


Bose- Einstein 


(16.71 Mev) 3 (Z/W) - 16 


-5.50 


Fermi-Dirac 


(15.65 Mev) 3 (Z/10) aie 


-2.74 


Maxwell-Boltzmann 


(14.73 Mev) 3 (Z/10) ' 16 


-1.83 



TABLE I: Values obtained of the surface density pod = ru n pu n (0) for Fermions and Bosons decoupling ultrarelativistically and 
for particles decoupling non-relativistically (Maxwell-Boltzmann statistics). [The exponent of Z originates in the primordial 
power n s /6 = 0.16]. The comparison of these theoretical results for poo = ru n pu n (0) with the observational value eq, (|1.2[) 
indicates that Z ~ 10 — 100 and therefore that the DM particle mass is in the keV range [see eq. (|2.25| l]. In any case, the 
agreement between the linear theory and the observations is remarkable. 



observed value eq. ljl.ljl . One representative solution is enough to obtain the value of the surface density but a 
more general treatment for non-spherically symmetrical solutions of the non-linear Boltzmann-Vlasov equation (and 
including also baryonic matter) would be necessary to prove the universality of ro /?(0). 

We can estimate the mass of the galaxies obtained in the linear approximation from eas. (|2.27)1 - ()2.28|) as 

4-n s 

( 10\ 6 4 — n 

M ~ rf in Plm (0) = 1.8 10 14 M f —J , — ^ ~ 0.506 . 

We obtain mass values in the upper range of the observations, as expected. 

Notice the scaling of the linear profile pu n (r) obtained here eq. (|2.28[) with the primordial spectral index n s : pu n (r) 
decreases as j» _1_n s/2 = r -i.482 j or r riin _ xhe value of this exponent is in agreement with the universal empirical 
behaviour recently put forward in ref.[6j r ~ 16±0A . For larger scales we would expect that the contribution from small 
k modes where nonlinear effects are dominant will give the customary r~ 3 tail. 

We find that the surface density computed from the linearized Boltzmann-Vlasov equation reproduces very well 
the observed value of the energy scale cq. (|1.2[) depending on the particle statistics. Nonlinear effects should improve 
the theoretical values of the surface density Pod — Hin pim{0) m Table I including the contributions from large scale 
modes. Notice from eq. (|2.29p that the theoretical poo has a mild dependence on Z, the only parameter here which 
is not known with precision. 

Anyhow, the agreement between the linear theory and the observations is already remarkable. The comparison 
of our theoretical values for p$o displayed in Table I and the observational value eq. (|1.2p indicates that Z ~ 10— 100. 
This implies that the DM particle mass is in the keV range eq. (|2.23p and (|2.25p . 

Notice that ru n in eq. (|2.27p decreases with Z as Z~3 7 while pu n {0) in eq. (|2.29p grows with Z as Z(™=+ 2 )/ 6 \nZ = 
Z 0A93 \nZ. Cusp profiles would thus correspond to the unphysical limit Z — > oo. The DM particle mass m grows as 
Z3 according to eq. (|2.23p (or as zi for DM particles decoupling non-relativistically HU). Cusped profiles are thus 
associated to heavy DM particles with a huge mass m well above the physical keV scale. For example, wimps at 
m = 100 GeV and T d = 5 GeV [27| would require Z ~ 10 24 Such heavy DM candidates are in conflict with the 
observed surface density eq. (jl.2|) . since eq. (|1.2|) indicates that Z ~ 10 — 100 and therefore that the DM particle mass 
is in the keV range [see eq. (l2.25p ]. Therefore, our present results for pojj = ru n pu n (0) provide further evidence for 
the mass scale of the DM particle being in the keV scale. 
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